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SUMMARY 


A procedure for determining the scattered pressure field resulting from a 
monochromatic harmonic wave that is incident upon a layered energy absorbing 
structure is treated. The situation where the structure is modeled with finite 
elements and the surrounding acoustic medium (water or air) is represented with 
either acoustic finite elements — or — some type of boundary integral 
formulation, is considered. Finite element modeling problems arise when the 
construction of the structure, at *the fluid-structure interface, are 
nonhomogeneous and in particular when the inhomogeneities are small relative to 
the acoustic wave length. An approximate procedure is presented for replacing 
the detailed microscopic representation of the layered surface configuration 
with an equivalent simple surface impedance finite element, which is especially 
designed to work only at limited frequencies. An example problem is presented 
using NASTRAN, however the procedure is general enough to adapt to practically 
any finite element code having a steady state option. 


INTRODUCTION 

This paper addresses the topic of solving acoustic-structure interaction 
problems involving a configuration having some sort of soft layered 
viscoelastic material, applied to a stiffer main body. Typically, acoustic 
sound energy impinges upon the treated structure (e.g. Figure la) and it is of 
interest to determine the reflected acoustic pressure. Specifically, the 
finite element method of solution is considered for the representation of the 
acoustically surface treated structure, and either acoustic finite elements 
(ref. 6,7)--or--some type of boundary integral method (ref. 8) that does not 
directly involve modeling the fluid is considered for the fluid. For either 
type of fluid representation there still remains the difficult problem of 
representing the energy absorbing properties of the viscoelastic outer layer in 
cases where the microscopic details are too complicated to represent with 
finite elements in the practical case where a large region of the structure is 
to be analyzed. On the other hand if only a small patch of the structure were 
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FIGURE • 1 GENERATION OF MACRO DYNAMIC STIFFNESS (R) 
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considered of, say, area A p (see Figure la), it would be well within the range 
of practical finite element analysis. 

The idea pursued here is to use this microscopic patch, in some sense, to 
generate a super finite element with many less degrees-of-freedom as the 
detailed microscopic patch, yet produces nearly the same reflected pressure 
field had the main body been entirely been modeled with the detailed finite 
element representation. Essentially the finer details of the soft 
nonhomogeneous layer are smeared out and represented by a sequence of simplier 
macroscopic elements. Conceptional ly, the Figure la configuration is replaced 
by the simplier finite element model shown in Figure 2, where the concept is 
generalized to apply to a curved surface. In Figure 2, the details of the 
stiff backing layer are not shown and are left as a choice to the model maker 
whether to represent the backing structure with say plate elements or solid 
brick elements. A choice of which acoustic fluid representation (i.e. finite 
element or boundary integral method) must also be made. The remainder of this 
paper focuses on the problem of defining the properties of the microscopic 
viscoelastic nonhomogeneous layer by specifically prescribing an equivalent 
macroscopic lumped parameter element of the type shown in Figures 3 and 4. 

Obtaining an approximate simplier model to represent a more complex 
nonhomogeneous viscoelastic layer has been considered by others, refs. (1-5). 
However in these earlier works, particularly ref. (1 and 5), the emphasis was 
concerned with obtaining a simple lumped parameter representation of the layer 
directly in terms of the identifiable physical parameters of the layer (e.g. 
modulus of elasticity and thickness) so that the physics of the layer's 
acoustical performance could be explained. In contrast, for the work presented 
here, we use the lumped parameter only as a means towards fitting the actual 
microscopic surface impedance, and the meaning of the lumped parameters need 
not be related to any specific physical properties of the nonhomogeneous layer. 
Further, in none of the referenced works (1-5), has an attempt been made to 
apply the results of the work to some sort of finite element scheme such as in 
Figure 2. 


OUTLINE OF THE METHOD 


Before proceeding to the detailed development, a brief outline of the 
procedure is helpful. The first step is to cut out a representati ve patch, of 
surface area, A p , from the overall configuration, such as the one shown in 
Figure la. Next a detailed finite element model of the patch is constructed, 
like shown in Figure lb, where from this finite element model of the soft 
nonhomogeneous layer alone, we can obtain the smeared out macroscopic dynamic 
stiffness. Next, through a curve fitting process, an equivalent simple lumped 
parameter finite element model is designed that has the same surface impedance 
as the finite element patch at only one frequency (for the Figure 3 model) or 
at only two frequencies (for the Figure 4 model). Next the simple lumped 
parameter elements are distributed over the whole surface of the structure 
(e.g. Figure 2), where the same lumped parameters are used for all surface 
nodes, except for an alteration accounting for the fact that the surface area, 
A n , around each node might be different if a variable mesh is used. The fact 
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FIGURE - 2 SINGLE ELEMENT SMEARED MACROSCOPIC REPRESENTATION OF SOFT LAYER 



that the direction of the normal to the surface changes on a curved surface is 
treated automatically when the bar version (NASTRAN CONROD) of the lumped 
parameter model (Figure 3C or Figure 4C) is used. Finally, the fluid is 
connected to the fluid side terminals of the lumped parameter elements (either 
acoustic finite elements— or— the boundary integral interaction matrix), and 
the scattering problem is solved in the usual manner for scattering from a 
homogeneous elastic body (e.g. using NASTRAN rigid format 8 with the acoustic 
elements of ref. 6 or 7 as appropriate and absorbing radiation boundary of ref. 
9). 


MACROSCOPIC IMPEDANCE DETERMINATION 


The first step is to obtain the macroscopic surface impedance of a small 
patch of the actual nonhomogeneous viscoelastic layer. The methodology given 
here is general enough to treat a variety of inhomogeneities ranging from the 
tubes of ref. (3) to the imbedded voids used in refs. (2) and (5). Therefore 
in what follows, we are not specific about the layer, wherein the only 
requirement is that the patch size, A p , is just large enough to pick up a 
typical repeated pattern. The Figure la, shows the sample patch as a WxD 
rectangular one, however, the shape should be whatever is convenient to 
represent the repeated pattern. Figure lb shows a generic finite element 
representation. It is important to point out that the details of the finite 
element model must be fine enough to properly represent the complex spatial 
response existing within the layer. It is further assumed that the Figure la 
viscoelastic layer is infinite in extent and that plane of symmetry type 
boundary conditions can be applied to all four lateral faces (of areas WxT and 
DxT). This boundary condition is represented by the zero lateral constraint 
indicated by the rollers shown in Figure lb. 


The Figure lb model is designed to respond to normal pressure that is 
approximately uniform over the patch surface. To obtain the smeared effects of 
the inhomogeneities, a weightless rigid "piston like" member is attached to the 
left and right surfaces of the layer (hence the designation L and R to denote 
the left and right surface nodes and all other pertinent surface related 
items). The attachment of the pistons to the layer is a rigid connection in 
the normal direction of motion, however the lateral direction attachment 
depends on the specifics of the application. The usual case is also a rigid 
attachment in a direction transverse to the motion direction on the left side 
(attachment to stiff backing structure), where as the attachment can allow 
transverse slip on the right side (attachment to the fluid). Double nodes are 
required when such slip is allowed. 


The relationship relating the left and right piston faces and displacements 
are given by 
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where the forces, F L and F R , on the left and right piston vary harmonically 
according to 

F L 3 F l exp(itut), F r = F r exp(iwt) (2) 


and the piston displacements Ul » Ur according to 

U L = 0 l exp(iojt) Ur= URexp(iut) (3) 

The 2x2 dynamic stiffness, [K]» which is the coefficient of {U^ in equation 
(1), can be computed from the fundamental matrices [KJ, [Mj, and I.CJ of the 
figure (lb) microscopic finite element model by decomposing the response into 
boundary nodes (B subscripts) and interior nodes (I subscripts). Thus 

.V : _ _ Kb j 

*IB » K I I. 

LKJ t 


(4) 



can be used to solve for {F}gin terms of (U]g, (note that {F}j ^Olsince there 
are no forces on the internal nodes), where the elements of [K] are of the 
form 


CKJ t = (-CMJw + [Kj + iu CCJ) 

It follows that after eliminating (Ujjfrom (4), we obtain 


(F }, 


= LK nn - K nT Kjj KjgJ{U} c 


(5) 


'B " U, 'BB "BI "II "IB JtvJJ B 
and comparing equations (1) and (5), it is evident that CKJ can be computed from 


LKJ - [K bb - K BI Kj] Kj B J (6) 

Therefore, one approach to obtaining the [KJ representation of the microscopic 
layer at some desired frequency, say u>. , would be to evaluate [KJ by 
constructing the right hand side of equation (6) with a OMAP sequence of 
instructions. This however_ is an unwieldly process which involves both 
partitioning_the microscopic [KJj matrix (e.g. ref. 10) and forming the 
inverse of [KJjj, which could be a costly process. 

An alternate process for generating [Kj is to compute it directly from two 
finite element runs, involving the Figure lb patch model. For a single _ 
frequency w,, we need to solve for the four complex constants comprising [Kj 
in equation 1 (1). These four constants are generated from the following two 
finite elements runs which are referred to as computer "run-a" and "run-b" 


• "run-a"; set 1.0, uf» 0.0, compute NOTE: 

L L Superscripts a and b 

U? , F* from NASTRAN "run-a" refer to "run-a"and 

L R "run-b" 
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"mn-b"; set F J? = 1.0, 0 b = 0.0, compute 

R L 

0 b , f£ from NASTRAN "run-b" 

Where the complex forces, and F? are read from the forces of constraint 
printout as activated by the presence of SPCFORCES = ALL card in the case 
control in the case of NASTRAN. The 0® and Ur complex displacement 
amplitudes are part of the normal finite element output as activated by the 
DISPLACEMENT = ALL card. Inserting the results of these two runs directly into 
equation (1) gives four equations and four unknowns to determine the four 
complex stiffness entries of the 2x2 [i<] matrix, namely 
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where the matrix [i<] is symmetric. The response and constraint forces will in 
general be complex when the microscopic finite element model of Figure lb 
consist of materials having viscoelastic dissipation (this effect is activated 
by the GE parameter on a NASTRAN MAT1 .card). For each frequency, to , there 
are six numbers K [l » K li_ , ... Krr that define the smeared macroscopic 
stiffness. Upon repeating the process described to generate equation (7) over 
a sweep of frequencies, a frequency description of the smeared macroscopic 
stiffness can be obtained like the one illustrates in Figure 6. Using the 
frequency sweep option in NASTRAN, the data used to generate the Figure 6 
example was generated with just two computer runs, by frequency sweeping the 
two "run-a", "run-b" cases described earlier. 


LUMPED PARAMETER ELEMENTS 

The three complex terms in equation (7) define the equivalent dynamic 
stiffness for a patch of area A normal to the surface, in a global coordinate 
system with one coordinate axis cilso in line with the normal to the surface. 
There now remains the task of implementing equation (7) as an element in the 
finite element code, so that the as yet undefined layer elements shown in 
Figure 2 can be implemented. Two basic lumped parameter models are developed 
herein. The first one is a single frequency model, that is designed to 
represent equation (7) only at one selected frequency. The second model is 
similar, except it is designed to represent equation (7) at two selected 
frequencies. Further, the second model automatically interpolates between the 
two selected frequency and therefore attempts to represent the nonhomogeneous 
layer response over a limited hand of frequencies. The area of the patch, A p , 
used to generate the macroscopic stiffness of equation (7) is in general, 
different from the area factor, Ap , used to convert distributed loads into 
concentrated nodal forces (see Figures 2 for A n ). Consequently, the 
generalized stiffness per unit patch area (l/A p )[Kj is a more fundamental 
quantity to work with. The dynamic stiffness to use at a typical node 
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connection shown in Figure 2 then would be 

LKJ = (A n /A p )CKJ 

Implementation of LKJ in the finite element code is awkward, although 
straight forward, when the element is not in alignment with one of the axis of 
the global coordinate system. In this case, equation (8) is further modified 
with a coordinate transformation 

I XT = CTJ T LK]LTJ O) 

where LTJ is the usual local to global coordinate transformation used in rod 
finite elements (ref. II). 

The implementation of equation (9) directly might require some sort of 
preprocessor, to generate the Lk’J' matrix in conjunction with the coordinate 
transformations. A simplier alternate procedure is also given, that redefines 
the local LKJ matrix in terms of rods (tension members such as CONROD and CVISC 
elements in NASTRAN), and the existing machinery in the finite element code 
will automatically take care of any transformations. Still another way to treat 
the coordinate transformation is to let each node have its own coordinate 
system, where one of the axis is aligned along the normal to the surface. This 
can be invoked in NASTRAN via a C0R01R card. In this manner, NASTRAN would 
automatically take care of the transformations, however it would require every 
node to have its own coordinate system. 


Single Frequency Lumped Parameter Element 

Here we are concerned with implementing equation (7), to each surface node 
(a.g. as in Figure 2), but we only require that the analysis be done at one 
frequency, uj. We give three alternate methods of installing the appropriate 
nodal dynamic stiffness. 


a) Direct Matrix Entry Version 

For programs such as NASTRAN that accept direct entries to the assembled 
global dynamic stiffness matrix (LKjJ, in NASTRAN) the components of equation 
(9) can be entered directly with DMIG cards. Since NASTRAN accepts complex 
entries, there is no problem with inserting both the real and imaginary parts 
of the LKJ' matrix into NASTRAN. 

For finite element codes other than NASTRAN that do not accept a complex 
stiffjjess entry directly.it can be done indirectly by inserting the real part 
of LKJ' through the usual structural stiffness entry, LKJ, and the imaginary 
part of LKJ' through the damping matrix, LCJ, after dividing img.pt(LKj') by 


b) Spring-Damper Entry Version 

The matrix LKJ, (before coordinate transformation) can be built from simple 
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lumped parameter springs and damper as illustrated in Figure 3b. The idea is 
to assign values to the springs and dampers, so that the assembled [KJ matrix 
is formed. By inspection of Figure 3b, the local element dynamic stiffness is 


LKJ = 


(•W *'“«* +c t> 


-K - i co C 
m m 


SYMM 


(K + KJ + iw (C m + CJ 
' r m m r 


(10) 


Upon equating equations (8) and (10) in conjunction with equation (7), we 
arrive at the following constants for the single frequency lumped parameter 
model : 


C m = 

K l = (\l + W* 

where A = A n /Ap, is the nodal area to patch area ratio. 

In the NASTRAN program, these lumped parameters are entered via CELAS2 and 
CDAMP2 cards. When the elements are not in line with the global axis, 
coordinate transformations are involved, i.e. involving equation (9)— or-- 
employing different reference frames for each element using C0R01R cards. 


= 0<LL + k lr)A/^ 

(^rr + y* 

= <+ < R ) A-/03 


c.) Three Rod Version 


The three rod element (made from C0NR0D elements in NASTRAN) shown in 
Figure 3b is an alternate procedure to install the desired equation (9) dynamic 
stiffness. The rods, properly sized, will generate the same dynamic stiffness 
matrices as either of the two mentioned ones, but has the advantage that the 
coordinate transformations, should they be needed, are taken into account 
automatically in the finite element program. The same dynamic stiffness 
realized by the Figure 3b model can be achieved with the Figure 3c three 
massless rod model by using the following rod properties shown in Table 1. 


TABLE I Single Frequency Rod Properties 



Mass 

Density p 

Young's 
Modulus E 

Poisson's 
Ratio v 

Loss 
Factor n 

Cross Gee, 
Area A 

Rod' 

Length L 

Left rod 

0 


0 

u£ £ /K £ 

T 

T 

Middl e rod 

0 

K m 

o 


T 

T 

Right rod 

0 

K r 

0 

V K r 

T 

T 


The stiffness of a rod is AE/L, and the damping matrix is n AE/L, thus it can 
be seen how the Figure 3b and 3c models result in equivalent dynamic stiffness 
matrices. The 6 Constants K&,. . .C r are still computed using Equation 11. 
The material properties p , E,v , n» A appear on the CONROO and MAT1 cards in 
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a) DIRECT MATRIX ENTRY VERSION 



b) SPRING - DAMPER VERSION 


c) THREE MASSLESS RODS VERSION 


FIGURE • 3 SINGLE FREQUENCY IMPEDENCE ELEMENT 



NASTRAN, however length T is computed by the program as the distance between 
terminal nodes R, and L in Figure 2. The user must be sure to enter a T in the 
area slot that is consistent with this length! Also care must be taken to 
ground out one end of the left and right rod as shown schematically in Figure 
3c. See Table 2 comments for entering a negative modulus, E. 


Two-Frequency Lumped Parameter Element 

The object of this macroscopic element is to have a single element 
represent the dynamic stiffness of the patch over a limited range of 
frequencies, without having to change the parameter constants. This would have 
use as a convenience feature when performing a frequency sweep 
analysis--or--may be of use in a transient type solution, where the model is 
expected to be accurate over a narrow band of frequencies. 

a.) Mass-Spring-Damper Version 


The element local dynamic stiffness matrix for the model shown in Figure 4b 
is given by the relation 


LKJ 


H“<V C Z> +K £ tK m) ! 
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(K r + K m ) + H 6| ) + 
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( 12 ) 


The structure of equation (12) permits one to vary both the real and 
imaginary parts of the three main entries of the [KJ matrix by changing the 
frequency oj . Note that had a lumped mass been used for the center mass 

element, the real part of the K gentry would not depend on frequency. A set 

of twelve equations and twelve unknowns can be set up to solve for twelve 
unknown lumped parameters K*, ....M m appearing in equation (12). The system 
of equations is obtained by equating corresponding real and imaginary parts of 
equation (12) and (8) at the first desired frequency, u, and again equating 

equations (12) and (8) evaluated at the second desired frequency,^ . Let the 

entries of equation (7) with a subscript 1 refer to the "run-a", "run,-b" pay 
of stiffness generating runs described earlier at frequency (i.e. K Tq ) 

and similarly let entries of equation (7) with subscripts 2 refer to the 
"run-a", "run-b" pair run at frequency^ i .e. kT^ >***^r 2 . • Solving these 
twelve weakly coupled equations, the following twelve results 
the unknown lumped parameters, 


are obtained for 
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where 
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It is noted that Bj and are not unknowns, but rather free parameters that 
can be selected by the user. However, when using the six rod version discussed 
later, one must use the 6 n , 62 factors corresponding to the consistent mass 
matrix being employed by the finite element program being used. The 
implementation of the model shown in Figure (4b) cannot be achieved, due to the 
fact we need a complex spring, iK|j, which is not possible to enter with say an 
CELAS2 card in NASTRAN. Instead, r it can be realized either with the direct 
matrix entry option or with a six rod element shown in Figure 4c as explained 
next. 


b) Direct Matrix Entry Version 

Once the basic lumped parameters have been computed, via equation (13), the 
element matrix in equation (12) can be entered directly with DMIG cards for the 
assembled global NASTRAN mat ri x, ..B^ andt^] matrices. Again letting the 
notation L, R denote the left right terminal node notation, we have 

- «.♦ K ), iOtf * r) 


dd LL ~ " - 

' K dd]/ -V -< 

'“dd^- ( V i(K r * K? 


(enter as complex on DMIG card) 


' B dd) = C m * C «. 

LL 


(14) 


tBdd L' " c - 


'"dd* - C «,* C r 

RR 


[M ddj, =a ^ a » ei 

LL 1 


I M ddl = \ S 2 

LR 


[M. ,] = M + M 8i 
ddi„ r m 
RR 

As in the previous single frequency model, step to account for the 
coordinate transformation must be taken where applicable. The next 6 -rod model 
will handle the coordinate transformations automatically. 

c) Six Rod Element Version 


A constant mass viscoelastic rod element with viscoelastic damping, n , has 
a dynamic stiffness matrix of the form: 




od 


9 9 

(AE/ Z -mu) 8, )+inAE/Jt I -(AE/S.+nno So) -iHAE/H 

1 I C 


SYMM 


! (AE-mo) 2 B 1 ) +iHAE/Jl 


\ ilJ-i U.UJJ b' 1 . 

r* 1 


(15) 


where A is the rod cross section area; m=A P& the total rod mass; £ is the rod 
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SHADED - CONSISTANT MASS VISCOELASTIC ROD 
UNSHADED — VISCOUS DAMPER ROD 


a) DIRECT MATRIX ENTRY VERSION 


b) MASS -SPRING -DAMPER VERSION 


c) SIX - RODS VERSION 


FIGURE - 4 TWO-FREQUENCY IMPEDENCE ELEMENT 



length, E the modulus of elasticity and Sj , 8 2 are the mass distribution 
factors. It should be noted that the 6 factors may be different versions of 
NASTRAN (see note below Table 2). 

Upon comparing equation (14) to equation (15), a set of equivalent rods can be 
defined to correspond exactly to the Figure 4b configuration. Table 2 below 
give the equivalent constants 


TABLE 2 Two Frequency Rod Properties 



Lumped Roc 

Mass Young* s Poisson* s Loss Damping Cross Length 

density Modulus Ratio Factor Constant Sectional (Indirect 

p E v n C Area A input) L 

left viscoelastic rod 
left damper rod 
mid viscoelastic rod 
mid damper rod 
right viscoelastic rod 
right damper rod 

h / l.T~ ' ^ 0 K^/K £ - T T 

M /T 2 K 0 K C /K — T T 

mm mm 

^m ^ 

M /(/,D K 0 K C /K T I 

r 1 r r r 

L r 


• For COSMIC NATRAN 3, * 1/3 3 ? 3 1/6 

NOTE: L 

• For MSC NASTRAN 3 x = 5/12 3 2 = 1/12 

In both COSMIC or MSC NASTRAN, the viscoelastic rods are implemented with 
C0NR00 elements (with GE as the loss factor on a MAT1 card) and the damper rods 
are implemented with CVISC elements where the lumped parameter damping constant 
is directly entered on the corresponding PVISC property card. The twelve basic 
constants K m ...M r given by equation (13) are the raw data that make up the 5 
rod model. The rod length is not direct input to the 6 rods, but is rather 
computed automatically by the finite element program based on the length 
between "L" node and “R“ node. It is important that the, T, length factor us* 1 , 
for the dummy rod density and dummy rod area in Table 2 must be totally 
consistent with the distance between nodes L and R. Also care must be taken tc 
ground out one end of the left OJMROD and CVISC elements and also ne end of 
the right CONRQD and CVISC slenents. Should E be negative, leave ic blank and 
inscead use a dummy pos ic ive shear modulus, G, and corresponding negacive 
"oisson's ratio chat corresponds co the desired negacive E. For example. 

■ et G = | E | and sec v = -1-5. 
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DEMONSTRATION EXAMPLE PROBLEM 

A soft inhomogeneous viscoelastic layer is attached to a stiff backing 
plate as shown in Figure 5a. The submerged sheet is a generic configuration 
that does not correspond to any particular real application, and is include 
here only to give some idea of the accuracy of the solution methodology in a 
simple application. Further, the details of the model are too lengthy to give 
here, and are beyond the intent of this demonstration problem. The example 
problem is kept simple by assuming the inhomogeneities are planar, thus 
permitting us to use two dimensional elements having a unit depth for the j 
variable of Figure la. The inhomogeneities were simulated by simply assigning 
different material constants (Young's modulus, mass density and dissipation 
loss factor) to soft layer finite elements in Figure 5b in the region marxed 
soft layer. The inhomogeneities were distributed mainly in nonlayered 
patterns . 

In the sample problem, the two frequency element of Figure 4 was used. T he 
TxW region in Figure 5b, was used to generate the macroscopic dynamic stiffness 
versus frequencies using the "run-a“, "run-b" procedure described earlier in 
the paper, in conjunction with the stiffness generator procedure illustrated in 
Figure lb. The results of these 5x2 runs are shown in Figure 6, where the 
frequencies, w, are normalized (divided) by (C a /W) where C, is the incident 
side acoustic medium wave speed. Further, the dynamic stiffnesses are 
normalized (divided) by the patch area, Ap = WO . 

We arbitrarily selected the Figure 4 two frequency model to work at 
normalized frequencies of .733 and .837, as shown by the triangle markers in 
Figure 6. We purposely did not make the spread greater, so that some idea of 
how well the model would work in between the forced points could be examined.^ 
For example, the dashed line shown connecting the two forced points in Figure 5 
(Krr plot) represents the dynamic stiffness of the Krr component of the Figure 
4 two-frequency model. Since we have not forced the dynamic stiffness to be 
exactly equal to the patch model, some error can be expected in the response 
results if the two frequency model is used at the in between non forced 
frequency points. 

The direct stiffness, mass, and damping via DMIG cards option (i.e. with 
equations 14) was used to implement the Figure 4 two-frequency model. The 
6-bars model was not tried, because the CVISC viscous damping elements do not 
always assemble properly on the DEC VAX SERIES MOOEL 11/780, April, 1984 
release of NASTRAN. Hopefully these elements will be fixed in future releases. 

The Figure 5b finite model was subject to an incident harmonic wave of 
strength P at the two normalized frequencies of .733 and .837 corresponding to 
the forced 3 frequencies of the Figure 4 .two-frequency model. The finite element 
model was first run with the full microscopic finite element representation for 
the soft layer (e.g. the soft inhomogeneous layer is represented with 398 
quadrilateral and triangular elements )— and — secondly was rerun with one single 
Figure 4 type element replacing the original 398 finite elements. The 
reflected pressure at the far end of the finite element model was used to 
compare results. The results are shown in Table 3. 
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a) SUBMERGED FLAT PLATE 



• 856 ELEMENTS (CQDMEM, (TRMEM, CBAR)) (ENTIRE MODEL) 

• 588 GRID POINTS 

W- WIDTH OF REPRESENTATIVE MICROSCOPIC PATCH OF SOFT LAYER 
b) FINITE ELEMENT MODEL 

FIGURE -5 SAMPLE NON HOMOGENOUS MICROSCOPIC FINITE ELEMENT MODEL 




WAp DYNAMIC STIFFNESS PER UNIT 
PATCH AREA 


INJ 

00 


(a) 


• -PATCH FINITE ELEMENT RUNS 
A-FORCED FREQUENCY FIT 


750,000 k 


500,000 b 




K 


i 

RR 




250,000 h 


'RR 


x 


TWO FREQ. MODEL 
A P (TYPICAL) IN BETWEEN 
FORCED POINTS 




a "77^X9 


Q_ 

1.0 


V 


NONDIMENSION AL FREQUENCY, kW 




FIGURE - 6 MACROSCOPIC DYNAMIC STIFFNESS VS. NONDIMENSIONAL FREQUENCY 



TABLE 3 Comparison of Results (REFLECTED PRESSURE/ 



The method presented provides a simple means for representing complicated 
finite element models of inhomogeneous viscoelastic layers with simplier 
elements, that are designed to work at specific frequencies. Preliminary 
results for the two-frequency macroscopic element on the simple flat sheet 
showed good results for the reflected pressure when compared to the same 
analysis using a full blown microscopic finite element implementation. The 
results were not quite as good, 23.9% error, for the comparison involving the 
"in between frequency run". However, the macroscopic elements are not 
specifically designed to work at the intermediate frequencies, and therefore it 
would be a windfall situation had they worked there as well. For example if 
the spread between the two frequencies were made larger, where the second 
forced frequency point in Figure 6 had been say, 1.0 instead of the .337, a 
substantia^ misfit in the dynamic stiffness would be expected due to the curved 
shape of i<^ , particularly in the range 0.8<kW<0.9. If results are desired in 
this range, new lumped parameters should be recomputed and a new 2-frequency 
model be used that is valid for the desired frequencies. 

The next future checkout application of the macroscopic elements should be 
for a curved surface of the type shown in Figure 2. It must be remembered that 
the macroscopic elements are not intended to represent the physics of the 
actual dynamic system existing between the L, R terminal nodes but rather 
simply represent the impedance of the actual inhomogeneous soft layer at those 
two frequencies. 
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